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^ ■ Abstract 

The properties of the Isoscalar Giant Dipole Resonance (ISGDR) and its electromagnetic struc- 
ture are investigated within a semiclassical nuclear Fermi-Fluid dynamical approach. Microscopical 
' calculations pointing to a ISGDR distribution splitted into two main broad structures is confirmed 

within the presented macroscopic approach by the occurence of a "low- lying" and a "high- lying" 

> , 

CN , state, that are nothing else than the first two overtones of the same resonance. Macroscopically 

' they are pictured as a combination of compressional and vortical nuclear flows. In the second part 
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of the paper the electromagnetic structure of the ISGDR, relevant for reactions with inelastically 
scattered electrons, and the relation between the vorticity and the toroidal dipole moment is an- 
alyzed. The relative strengths of the compresional and vortical collective currents is evaluated by 



^ ■ means of electron-scattering sum-rules. 
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I. INTRODUCTION 

Recently the College Station group \^ reported experimental results on the isoscalar El 
strengths in three proton magical nuclei (^°Zr, ^^^Sn and ^''^Pb) using inelastic scattering 
of 240 MeV a particles at small angles. The authors concluded that the isoscalar El 
strength distribution in each nucleus is shared mainly between two components, one located 
at low energy and and another one at higher energy. In a subsequent publication this group 
presented new data on the ISGDR [2]. For ^^^Sn, ^^^Sm and ^°^Pb the low-energy peak fall 
in the interval (1.71 — 1.92)7ja; whereas the high-energy peak lays between Shw and 3.2hw. 
The upper component covers approximately 3 times more of the energy-weighted sum rule 
compared to the lower component. Similar values for the two peaks for ^osp]^ g^j^g given in 
PI : 1.80^ and 3.25hw. Previously Morsch et al. 0] found for the high-lying component 
in 208p|j ^ centroid of 21.3 ± O.SMeV which corresponds to (3.15 ± 0.12)/ia;. Therefore, 
experimentally, the lower-energy component has a value very to the IVGDR centroid which 
lays around ~ 2^, whereas the higher-energy component is located in the same region as 
the electric octupole resonance, i.e. ~ ?>tvjj. 

On the theoretical side there have been numerous studies aiming to disclose the features 
of these exotic modes. The usually accepted macoscopical picture of the ISGDR is a "hy- 
drodynamical density oscillation" in which the volume of the nucleus remains constant and 
the proton-neutron fluid oscillates in phase back and forth through the nucleus in the form 
of a compression mode 0]. Microscopical calculations using strengths associated with 
the non-isotropic compression, namely the "dipole squeezing" operator D = 
are stressing the importance of the high-lying ISGDR (^ 3fia;) J?, ^, 1^. Th ey are also 
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predicting a rather fragmented peak at smaller energies (~ 2huj). 
bimodal structure, was obtained within the fully consistent relativistic Hartree-Fock plus 
RPA framework [ij] 

The first macroscopically based models of electric resonances, were primarly concerned 
to describe the ISGDR as density fluctuations of a nuclear liquid drop. In an apparently 
overlooked paper, available only in german and published a couple of years after the 
emergency of the incompressible fluid model of the IVGDR, for the first time the isoscalar 
dipole eigenfrequencies of a spherical nucleus jwere determined. In our days, the investi- 



gations of the College Station - Kyiw group 



16], pointed out that considering only com- 
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pressional components in the velocity field and neglecting the relaxation effects within the 
nuclear-fluid dynamics, leads to an overestimation of the energies of the 1~ resonances with 
respect to the experimental values. 

Other macroscopical approaches aimed at the description of the giant resonances (includ- 
ing the 1~, T = state) by allowing for vortical components along with or without the 
compressional(incompressional) of the velocity field. Deriving conservation equations, such 
as continuity equation, and equations of motion such as the Navier-Stokes or Lame it was 
shown that the macroscopic velocity field admits also shear (transverse) components jl7|. 
In ref.^^, after some simplifications compared to [3], that we are going to dismiss in the 
present study, an isoscalar 1~ state of pure vortical character was derived. Substracting the 
center-of-mass motion the associated velocity field reads 

V°'-(r, ^, 0) = V X V X r (^r^ - ^r{r^)^ Y^^{9, 0) (1) 

Most important, in this study, for the first time a connection to the toroidal class of elec- 
Tomagnetic multipole moments was done introduced earlier by Dubovik and Cheshkov 
3]- The theoretical search for a vortex-like isoscalar dipole electric excitation associated 
to the toroidal dipole momeiit (" dipole torus mode") was continued in the nuclear-fiuid dy- 
namics frame [2O, 0]- Ref.j^^] substantiated the elastic character of this mode, since a 
nucleus without shearing properties cannot withstand transverse-like oscillations, and eval- 
uate for the first time the (e, e') form factors corresponding to the excitation of this isoscalar 
1^ resonance. The radial part of the transverse electric form-factor corresponding to the 
electro-excitation of the dipole torus mode was found to vary like jz{qr)/qr. 

Apart from the quest of toroidal nature of the ISGDR there are also other issues re- 
lated to the enhancement of these electromagnetic transitions for other types of collective 
electric excitations. The electric dipole spin waves were identified to have non vanis hing 
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magnetization-dependent part of the toroidal dipole operator in ref.|22]. In 
the fingerprints of the dipole toroidal moments in the electromagnetic properties of nuclear 
rotational states were examined for the first time in the literature. In '24] it was inferred 
that the strong deviations from the predictions of the adiabatic theory for the absolute val- 
ues of ii^l-transitions in the Coulomb excitation of ^^^Ra are related to the enhancement of 
toroidal transitions between the ground state and the lowest negative parity band. 

The revival of the interest on the role played by toroidal moments in the excitation of the 
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isoscalar dipole resonance was caused by a recent publication '25] that aimed to evaluate the 
El strength distribution in spherical nuclei where the RMF formalism + RPA calculations 
were previously unable to provide a satisfactory agreeement with the experimental data on 
the positions of the ISGDR resonances. Using the toroidal dipole operator corrected for the 
cm. motion instead of the squeezing operator, broad resonant peaks were assigned in the 
low- (< 2;^, and Mgh-ene.gy (> 3«.) ~l the strength d.tHbut.on fo. -»Pb. 

Commenting on the results reported in [25||, the short-note [26|| stressed the fact that 



the low-lying vortical mode, as inferred from macroscopic calculations 
ent from the so-called "pigmy resonance" which la ys i n t 
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le vicinity of Ihu, and that the 
centroids provided by the nuclear-fluid approaches |2Ct |2l| are still providing a qualitative 
good agreement with the {a, a') scattering data. In this respect the merit of is that it 
offers for the first time in the literature a microscopical calculation of the toroidal content 
of ISGDR states and confirms the connection already established by macroscopic models 
between this electromagnetic characteristic and vorticity. Moreover, and this will be an 
important point in our present work, although not explicitly stated in the body of ref.j3] 
but rather inferred from Fig. 2, it indicated that also for the high-lying states there is a more 
or less important value of the toroidal strengths which implies that these excitations are not 
purely compressional but may contain also significative vortical admixtures. 

The nature of collective flows in a range of excitation energy below 20 MeV for 208p]^ ^g^g 
more closely approached in [2^, where calculations within the quasiparicle phonon model 
are pointing to strong vorticity below 2huj for the entire electric dipole response, not only 
the isoscalar one. 

I — I 

In a subsequent publication ^] the nuclear electric isoscalar dipole response was studied 
within the RPA formalism including in the dipole part of the separable interaction simul- 
taneously dipole-dipole {Fx^ = rYi^j) and compressional dipole (squeezing) {Fx^ = r^Yi^) 
fields. The authors infered that the isoscalar dipole response is less sensitive to these type of 
interactions and more to the single particle structure. The strength function for what they 
called the squeezed dipole mode, displays two main broad and fragmented peaks : a low- 
lying, which depends on the coupling constants of the Nilsson potential and has a centroid 
ranging between 8 and 11 MeV and a high- lying one with centroids ranging in the interval 
21-23 MeV. 

The present paper aims to a description of the properties of the ISGDR within a macro- 
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scopic model based on the nuclear Fermi-Fluid picture, where along with density fluctations, 
transverse components of the collective velocity field are included. The main purpose is to 
determine the share of compression and vorticity fiows in the states building the isoscalar 
dipole electric response. The electromagnetic structure of the ISGDR is studied for arbi- 
trary momentum transfer by using the multipole parametrization of charges and currents 
of ref.[3]. The relation between the transition electric dipole moments and the transition 
vorticity is discussed. Finally we introduce inelastic electron scatterig sum-rules in order to 
assess the strengths of the states building the ISGDR with varying momentum transfer. 



II. NUCLEAR-FLUID DYNAMIC APPROACH 



A macroscopic approach which goes partially on the lines already developed in a previous 
publication [21] is adopted. However a few strong amendments are performed : 

• full /c-content in the radial part of the collective field 

• a compressional elastic constant different from the shear elastic constant fALame 7^ 



AiLame) that draw it nearer to other Fermi-Fluid dynamical approaches 



|l6, 



17 



The procedure consists in taking moments of the Boltzmann equation, i.e. to integrate 
it in the momentum space with weights 1, pi, PiPj, etc.. The first two moments provide the 
continuity and the Navier-Stokes equation in view of their similar form to the well-known 
equations known from Hydrodynamics. 

dp{r,t) 



dt 



V-(p(r,t)7i) = (2) 



-(p(r,tK) + 5^-^ = (3) 



dt ^ dxj 

j=i J 

where p is the mass density, which is supposed to be of the sharp-edge type in the present 
work, u is the collective field (which vanish in the ground state), whereas Pij are stress 
tensor components. To these equations the linearization procedure is applied 

p(r, t) = po + ^p{r, t), u = 6s, P^j = p5ij + 5aij (4) 
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In the second equation of (0)) the displacement field, Ss, was introduced. The stress tensor 
is splitted in a diagonal part (normal pressure) 

p^-i-A-poI' (6) 

and a non-diagonal one (associated to the shear) 

In the above two formulas K is the incompressibility coefficient of nuclear matter and ei? is 
the Fermi energy. The scalar function 

D = V ■ 5s 

describes the compressibility of the displacement field 6s and Sij are the components of the 
dyadic strain tensor j^^, 

£ = 1 (Vs + sV) 
After linearizing the equations of motion we get 

Sp = poV (7) 

PqSs = (ALame + 2/iLamc) V( VP) - 2/iLameV X U: (8) 

where, like in Hydrodynamics (0|, p. 115), by uj we denote the vorticity vector which is 
proportional to the curl of the displacement field 

(jj = -V X OS 

The equation of motion (jH)) is identical to the Lame equation ([s^, p. 60 and 94) known 
in the Mechanics of deformable continua, where the Lame elastic coefficients are provided 

n 

by the properties of the nuclear Fermi gas 

_ noK 4 _ 2 

^Lame — ^^^0^^; A^Lame — -^^O^F 

For the incompressibility coefficient K we use the fact that the excitation energy of the 
Isoscalar Giant Monopole Resonance (ISGMR), which is an isotropic volume oscillation, 
can be related to the compressibility of nuclei j^], i.e. 

£^isGMR = (-^7^) ^ ~ 82 ■ A-'/' 



In what follows a fundamental theorem of vector analysis is used p. 131) which states 
that every continous vector field V, which, together with its derivative falls to at large 
distances can be decomposed into a divergenceless part V±(V ■ V± = 0) and a curless part 
V|l(V X V|| = 0). 

Then eq.® separates in a equation for the compresibility {V) and another one for the 
vorticity (a;) which describes the degree of shear of the displacement field for the case of an 
elastic body. 

V = clAV, u) = c^Auj (9) 

where cl = ■>/ ALame + 2yULame/ Po and ct = a/ /iLamo/ Po 3.16 the propagation velocities of the 
longitudinal(compressional) and transversal (shear) elastic waves in nuclear matter [2^. 

Assuming an harmonic variation in time of the fluctuating parts of the density and the 
displacement field, i.e. 




5p{r, t) = p{r)e'''\ 5s(r, t) = 5s(r)e*^" 

the compresibility and the vorticity are found to satisfy the scalar and vector Helmholz 
equation respectively (HE) 



(10) 



corresponding to the wave- numbers T = Q/cl^t- In seismology the compresional wave is 
called the P wave and the transverse wave the S wave j^^l- The S in its turn has two com- 
ponents : the SH wave (known as the "poloidal" in Hydrodynamics or "transverse electric" 
in Electrodynamics) and the SV wave ("torsional" or "magnetic"). For a nucleus with a 
sharp edge one adopts a spherical geometry and the radial part is given by spherical Bessel 
functions whereas the angular part can be written in terms of spherical harmonic vectors. 



33l |. Next we disregard the torsional 
and consider only axial-symmetric 



Details can be found in the appendix or in the literature 
component which is related to magnetic excitations |84 | 
displacement fields (/i = 0). We then have for the longitudinal and poloidal components of 
the displacement field the expressions ()A5|1 and ()A7jl derived in the appendix. These ex- 
pressions have to be further corrected in order account for the center-of-mass motion. Like 
in a preceeding paper j^| the translational invariance of the collective velocity field results 
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from the condition that the center-of-mass Rc.m. is at rest 

/ drp{r,t)6s 



5R 



'CM. 







jdrp{r,t) 

Thus in the dipole case (A = 1) the longitudinal and transverse displacement fields are 



joikLr) 



10 V^' ' 



V2j2{kLr)Y 



(11) 



;i2) 



V2 (^MkTr) - ^^j,(kTRo)^ Y%i9,<P)-j2ikTr)YU0,- 



(13) 



The corresponding corrected expression for the density fiuctuation results by applying the 
continuity equation in (fTTj) followed by the substitution of the longitudinal diplacement field 
fll2|) . The integral relation between the corrected expression of the density fiuctuation and 
the longitudinal displacement field reads 

J dr rSp = po J dr V X {Ssl x r) (14) 

Eventually, we get for the density fluctuation an expression identical to the one derived in 

Sp = apo (^Ji{kLr)Q{Ro - r) - l-j,(kLRo)S{Ro - r)^ Yio{0, 0) (15) 

In order to derive the longitudinal and transverse wave-numbers k^ and kx and the 
constants a and b, multiplying the displacements fields, boundary conditions for the force 
acting on the free surface of the nucleus have to be imposed. This force is obtained by 
projecting the dyadic stress tensor on the normal unit vector to the surface. 

F = P Br (16) 

Usually two types of boundary conditions are employed depending on what assumption 
has been made for the surface. Two kinds of bounding nuclear surfaces are distinguished 
for sharp-edge distributions : rigid surfaces on which no slip occurs (used in the liquid 
drop- model to determine the "surfon" eigenvalues or in the hydrodynamic model of giant 
resonances 3^ to determine the "gion" eigenvalues) and free surfaces on which no tangential 
stresses act. If we would assume a rigid surface then we would end up with a density 
fiuctuation, and thus also with a velocity field, containing admixtures of the center-of- 
mass(c.m.) motion. Actually the expression for the fiuctuation density (jTSJ, corrected 
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Overtone(fi) 


i,{")/f™-i\ 
kj^ [tm ) 




(MeV) h^ln/hujo rn 


— hi 

= Onl 


1 


2.05 


3.05 


11.56 


1.67 


1.94 


2 


3.93 


5.86 


22.19 


3.21 


-1.19 


3 


5.05 


7.53 


28.53 


4.12 


6.09 


4 


7.09 


10.57 


40.03 


5.78 


-4.03 



TABLE I: The first 4 overtones of provided by the eigenvalues of the boundary condition (|17j) and 
the vorticity /compressibility ratio for 208 pj-, 

for the cm. motion, is compatible with the assumption of a free surface. The boundar 



conditions on a free surface impose that the force fulfill the following two conditions 
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The above equations provides an infinity of eigenvibrations but for the study of giant reso- 
nances only the first few are relevant. The boundary condition ()17p allows also the determi- 
nation of the ratio &„/a„ which gives the admixture between the compressional (longitudinal) 
and vortical(transverse) field in a given state n 

U pLongitudinal 



^ pTransversal 

^rr r=_Ro 

In table |l] we list the first 4 roots(overtones) of the boundary condition plus the ratio of 
the transversal-to-longitudinal weights. 

The first three eigenfrequencies for the dipole density fiuctuations as derived in the pio- 
neering work of Woeste have values very close to the density-vorticity waves listed in 
the above table : 1.78 hw, 3.02 hu and 4.22 hu compared to 1.67 hw, 3.21 hw and 4.12 hw. 
Also the estimation for the T = 0, L = 1^ vortex mode from {1.7hw) is very close to the 
value derived in the present paper. In a previous paper dedicated to ISGDR j2l| together 
with collaborators we constrained the collective velocity field to admit purely vortical fiows 
and we considered the low-lying response of the Fermi liquid, whereas in the present ap- 
proach the full momentum content is taken into account. Therefore, although compressional 
components are occuring, the first overtone has a rather vortical character as can be easily 
inferred from the upper left panel of Fig^ and displays the typical Hill-vortex pattern as 
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FIG. 1: Flow lines corresponding to the first four overtones of the ISGDR in ^^^Pb. 



already mentioned in |18l . l2Ct l21l | : Nuclear matter flows around a vortex ring situated in the 
equatorial plane of the nucleus. Other works, based like 15] on the compressible and irro- 



tational nuclear liquid drop, e.g. [37[, fail to observe the first overtone. In such approaches 
the necessity to correct the density fluctuation for the cm. motion was not anticipated. For 
the second overtone the compressional and vortical flows have almost an equal importance. 



10 



which contradicts the entrenched picture of a compressional mode around Shu. However 
there is to the date no direct experimental indication on such a macroscopic property of this 
higher lying isoscalar dipole resonance, and therefore the result that we report on this mode 
should not be excluded from debate. To a certain extent the flow pattern of this mode (see 
upper right panel of FigC]) presents typical characteristics of a compressional mode : the 
concentration of nuclear matter flow inside the southern emisphere and the depletion inside 
the northern emisphere; in the same time an opposite behavior of the density fluctuation is 
manifested at the north and south poles. For the third and fourth overtones (see lower left 
and right panels of Fig^) the flow is predominantly vortical. Although difficult to disclose 
from the flow patterns, the vortex structure becomes more intricate in the sense that in- 
stead of one vortex ring as was the case for the first overtone, we deal with two rings (n = 3 
overtone) and respectively three rings (n = 4 overtone) of lower intensity. Note that two 
succesive rings have opposite rotational fiows. 

A quantitative way to asses the role of compressional and vortical fiows is to introduce, 
following the orientation averaged values of the collective velocity field divergence 



and vorticity 



These two quantities are displayed in Fig|^ In order to avoid the awkward effect on 
these two radial functions near the surface, which is due to the sharp-edge distribution, (see 
), the curves drawn in FigElwere computed by assuming a diffuse density distribution. 
Consequently an additional peak in both {T>) and (u:) occurs in the surface region. In 
what concerns the vorticity we remark for the first overtone that the vorticity attains a 
maximum at approximately which corresponds to the critical points of the Hill 

vortex, a fact already pointed out in [2^. Inside the nucleus the compressibility increases 
almost linearly with the radius and is less important than the vorticity, thus confirming the 
previous assignment of this collective state as dipole torus mode jla 211- When the 
overtone number increases the vortex with the largest strength migrates towards the center 
of the nucleus, and new ring-like vortices occur at larger radii. For the overtones with n> 2 



the compressibility developes maxima inside the nuclear sphere and it plays a dominant role 
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only for the n = 2 overtone (the claimed Shw "compression mode") in the vicinity of the 
nuclear surface. The mixed (compressional+vortical) character of the n = 2 state can be also 
inferred from the self-consistent RPA calculations with Skyrme type interaction performed 
in an old study [38] as well in the relativistic mean- field approach from f^. 

The procedure to quantize a continuum system, described by the equation of continuity 
and the equation of motion (jH| is to expand 6p and 6s in normal coordinates js^ '■ 

5p{r,t) = ^p„(r)a„(t), 6s{r,t) = ^s„(r)a„(t) (18) 

n n 

The above sums are running after the values of k^lkx) allowed by the boundary condition 
(I17|) . i.e. after the overtones n. For the expression of the kinetic energy in the newly 
introduced collective coordinates we have 

T = IpoJ rfr|5s(r,t)|2^1 (19) 
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FIG. 2: Compresibility and vorticity of the first four overtones of the ISGDR in Pb. 
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For the mass inertia parameter we derive the expression 
SmA 



Bn 



An 



I ^2 

I ' r). 



The stiffness coefficient associated to an oscillation of degree n is simply 

and the quantized form of the energy can be obtained by introducing the creation and 
annihilation dipole "gions" 

The collective Hamiltonian will then read 

H = Y,^nJdiL+^^ (21) 



III. ELECTROMAGNETIC STRUCTURE OF ISGDR 

A. Low-(7 limit of the form-factors multipole parametrization 

To disclose the structure of the ISGDR we adopt the multipolar parametrization of 



charges and currents according to jl9|. In this approach the classical electromagnetic multi- 
poles are expanded in the momentum transfer in reactions with photons, electrons or charged 
hadrons. Let us take first the charge multipole form-factor for the charge part of the ISGDR 
fluctuation density: 

t) = j drUqr)Y,,{^^)Sp,ir, t) ^ (^g,,(t) - 2(2aV 3) ^'^^^0 ^^^^ 

The first terms of the above g-expansion represents the transition charge dipole moment 

Qx^.it) = 6x,i6^,o J drr^Yx^ie, 0)5pp(r, t) = 0, (23) 

a quantity which vanishes for the ISGDR due to the constraint imposed on the cm. motion. 
For the IVGDR this will not be the case, since the djTiamical dipole moment arises naturally 
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(20) 



as a measure of the relative motion between the "negative" charge (neutron) distribution 
and the positive charge (proton) distribution. Instead the next term in the expansion (|^^. 
does not canceL The quantitiy 

S(t) = 5,,,5,,o [ drr'^'Y,,ie, 0)%(r, t) = 2p,,Rl J] a„(t)^^^^^|^ (24) 

represents the mean square radius of the dipole charge distribution. It provides infor- 
mations on the spatial extension of the ISGDR and it depends only on the longitudi- 
nal (compressional) part of the displacement field which are related via the continuity equa- 
tion (121) to the density fluctuations. 

n 

According to the charge-current multipole parametrization of , the electric transverse 
form factor, splits into the g = limit and a term containing the higher order content in q. 



7-^+^ / A -I- 1 / . \ 

TfM.t) = ^^X + iy/ 'l^ [Qx,m+q'Tl°;{q,t)) (25) 

The first term in the paranthesis is the time- derivative of the Coulomb multipole moment 
defined in eq. ()23|) The second term represents the toroidal form factor that reads in the 
low-g limit 



Tr(o.*) = -y^/*'- 



^^-i + 2A + 3V A + 1 ^"+1 



A 1 



j{r,t) (26) 

I cirr^+V^,-(Vx j(r,t)) (27) 



2z V A + 1 2A + 3 

In classical electrodynamics the transition toroidal multipole moment is associated to a 
poloidal flow on the wings of a toroidal solenoid (for details see the reviews 
study of electric collective states it can be related to the strength of the vorticity associated 
to a nuclear transition. Indeed, following ref. [4^, we introduce the transition multipoles of 
the curl of the current density (unconstrained by the charge-current conservation law) 

T,,{r)^{If\\ (Vx j(r',t))AA || /.) 

I — I 

In order to remove the charge-current conservation constraint, the authors of ^4(]j introduced 
the pure vorticity transition multipole 



uxx = Txy{r)-\/^np^ (28) 
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where p\ is the charge density multipole and VL the energy associated to the transition. If 
the quantity 

\x 

is defined according to J40!] as the strength of the vorticity and employing the definitions 
()24|) for the square of the dynamic dipole charge distribution and (jTTj) for the dynamic 
dipole toroidal moment we arrive at the expression relating the r.m.e. of these last two 
electromagnetic multipoles and the vorticity strength 

1 



{if = K \\Tr II ^^ = 0) = !^ 



^z/i + n„(/;= 1,, II qI II /i = 0) 



(29) 



From this last formula we see that since the toroidal dipole moment and the square radius 
of the dipole charge distribution are the leading terms in the g-expansion of the transverse 
electric ()25|1 and Coulomb (j22|) form factors for the ISGDR, the determination of these 
two electromagnetic multipoles at low q allows the determination of the vorticity content 
unconstrained by the charge-current conservation law. Before ending this section we give 
the classical expression of the transition toroidal dipole moment associated to the ISGDR. 
Since the current density reads in this case 

Z 

^0 

n=l 

we finally obtain 

Thus, the dipole toroidal moment depends only on the shear (vortical) part of the proton 
fiuid displacement field. 



j = 4''° ^ {^^L\r) + rjs^r\r)) (30) 



TnO,t) = ^Po.i?^E^"""W'-^pr^ (31) 



B. Electro-excitation form factors of ISGDR 



Inelastic electron scattering is an excellent tool to explore the nature of currents involved 
in the excitation of low-lying (rotational or vibrational) and high-lying (giant resonances) 
collective states 41j. A specific feature of this reaction is represented by the possibility to 
separate the longitudinal from the transverse response functions. This is of vital importance 
if one tries to disentangle the compresional from the vortical response in the excitation of a 
specific electric collective state. 
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C. (e,e') form factors 



Quantizing the expression of the fluctuation density (fT3j) and inserting it in the square of 
the reduced matrix element(r.m.e.) of the Coulomb multipole operator (j22|) leads to 

2 



{If 



II /. 



(32) 



where 



nO 



(n) 



— k 



{ny 



ki^'^UkrRoMqRo) - qUqRo)ji{kr'Ro) - J2{kr Ro)ji{qRo 



in), 



(33) 



and the amplitude of the canonical coordinate an is given by the transition r.m.e. of the 
normal coordinate 



ttnO = (-^/ = 1 II "n II h = 0) 



h 



1/2 



(34) 



It can be noticed from the content of that not only its low-g limit depend on the 
longitudinal components of the collective flow (Ss^), as substantietad by eq.(j21I), but also 
its higher order q content. Therefore, the Coulomb form-factor provides solely informations 
on the compresional waves associated to the ISGDR. 

In figOwe draw the dependence of (j32|) on the momentum transfer q for the first four 
overtones in ^''^Pb. The diffraction maxima, although shifted to higher values q for increasing 
overtone number, are not affected in absolute value. 

The expression of the squared r.m.e. of the transverse electric operator reduces to 

2 



a/ = i^ WT^iq) l|/. = o)p = 3 



(35) 



where 



kT ^ji{kTRo)jo{qRo) - qji{qRo)jo{kT'Ro 



(n) 



Ft (q) = rn^nOinO \ Rq 



3^{kPRo)3i{qRo) 
(36) 

Thus, contrary to the Coulomb multipoles, the transverse multipoles are providing informa- 
tions only on the vortical components of the ISGDR. The squared r.m.e. of is ploted 
in FigHJ The first diffraction bump dominates for the n =1 and 2 components of the IS- 
GDR, whereas for n=3 the second bump and for n = 4 the third bump are taking over. 
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FIG. 3: R.m.e. of the Coulomb multipole operator for the first 4 overtones of the ISGDR in ^'^^Pb. 
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FIG. 4: R.m.e. of the transverse multipole operator for the first 4 overtones of the ISGDR in 



According to TablelU were we listed the ratio 6„/a„(vorticity/compressibility), this feature, 
is a consequence of the vorticity enhancement, in this higher lying states. Instead the first 
Coulomb diffraction bump maintains its importance in the electro-excitation of higher-lying 
states because the role played by the density waves is less important for these states. 
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D. Sum rules for electron scattering 



It is known that for the isoscalar electric modes simulated by operators depending only 



on coordinates of the particles the energy-weighted sum rules can be determined mode 



in- 



dependently and depend only on the ground-state properties of the nucleus j5|, |6|, |37|, 142 1. 
Operators from this class are commuting with interactions not depending explicitely on the 
momenta of the particles. These scalar operators are simulating shape or density distortions 
corresponding to a given isoscalar multipolar resonance and for that reason in the litera- 
ture the macroscopic images associated to these excitations are always irrotational surface 
or bulk compressional oscillations. Instead they are inadequate to describe distortions of 
the nuclear current which are not constrained by the charge- current conservation law, i.e. 
excitations with vortical currents. A class of sum-rules coping with both kind of distortions, 
i.e. of charge density and current (unconstrained by the charge- current relations) is given 
by the longitudinal (L) and transverse (T) intrinsic energy weighted sum-rules (EWSR) 
at constant three-momentum transfer \q\, which are constructed by weighting the nuclear 
response functions j3| 

°° / fc2 2 \ 

i?^(q, E*) = J2 I I P(9) I 0) I' ^ - ^ - Enj (37) 



/ fc2 2 \ 

R'{q,E*) = J2\{n\J^{q)\0) \" ^ [E* ' ^ ' En) 



(38) 



with an appropriate power of the nuclear excitation energy E*. Summing over all excited 
states we obtain the intrinsic p-order EWSR depending on q 

/oo 
dE*'R'^{q, E*')E*^^ = Y.Ei{n \ p{q) \ 0) \' (39) 
n 

/oo 
dE*'R'{q,E*')E*'' = Y^E'Jn \ J^{q) \ 0) \' (40) 
n 

where E*' = E* - h^q^/2MA is the energy available for intrinsic excitations. In the above 
formula J±{q) denotes the transverse component of the current operator relative to the 
momentum transfer (Jj_(g) = J{q) — q{q ■ J)/q^). 

The longitudinal and transverse p— order energy strengths of each state can then be 
obtained as relative contributions to the corresponding sum-rules (j39p and (|4U|) 

El I (n I p{q) I 0) r El \ {n \ J^{q) \ 0) \' /m/{q) 
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Nucleus nii^ /rriQ^ ^Jm.'^^ jmi^ mi" /rriQ^ -y/ms^/mi"^ Kli KI2 (^^i)cxp (^^2)cxp 
(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) 
9°Zr 15.44 16.29 16.86 25.53 15.28 29.34 16.20±0.80 25.70±0.70^) 

"^Sn 14.19 14.99 15.51 23.68 14.04 26.96 14.38±0.25 25.50±0.60^) 

14.70±0.80 23.00±0.602) 

i^'^Sm 13.21 14.07 14.64 23.95 13.07 25.08 14.00±0.30 24.51±0.40^) 

208pb 11.68 12.44 12.96 21.19 11.56 22.19 13.26±0.30 22.20±0.30^) 

12.50±0.30 22.50±0.302) 

TABLE II: Ratios of ISGDR sum rules for q — > compared to the energies of the first two 
overtones for ^'^Zr, ^^^Sn, ^"^^Sm and ^'^'^Pb and with experimental latest experimental data from 
[3]^^ and P]^) for the low- and high-energy peaks . 



In FigsEJandiniwe represent the zeroth- and first- order strengths distributions for the lon- 
gitudinal and transverse responses. In the sums only the first three overtones were included 
and therefore the excitation energy is truncated at 30 MeV. We see that for very low mo- 
mentum transfer the L and T strength is mainly concentrated on the first overtone. When 
q increases the L and T strengths follow a different pattern. Whilst the L strength is frag- 
mented almost "democratically" over the three overtones, the T strengths are testifying a 
transition from a low-g regime where the first overtone dominates to a high-g regime where 
the third overtone overtakes the predominance. For both regimes the second overtone plays 
a very minor role. This fact can be explained in our view by the compressibility content of 
this mode which is "washed-out" in the transverse response function. In the longitudinal 
response the second overtone plays a more visible role. 

In the low-g limit the ratios mi /mo and A/wis/mi are expected to provide crude estimates 
for the mean excitation energies associated to the density or transverse density current 
distortions. In Table ITTl we list these mean excitation energies for a series of spherical nuclei. 
Comparing the obtained values to the eigenvalues of the first two overtones we notice that 
the ratio mi^ /m^^ gives the best estimate of the first overtone. For the second overtone the 
ratio yj'm^^jm^ gives a rather coarse approximate. The other two ratios are falling in the 
vicinity of the first overtone. 
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IV. SUMMARY AND CONCLUSION 



The lower and upper component of the ISGDR, as reported by the latest experimental 
measurements, are explained as the first two overtones of a spherical Fermi-Fluid system 
with a sharp free surface corresponding to a mixture of dipole compression and vorticity 
oscillations. 

The approach presented in this work was applied primary to the heavy spherical nucleus 
^°^Pb, because in this case the sharp-edge density distribution assumption is more acceptable 
as would be the case for lighter nuclei for which experimental data on ISGDR are available 
(^°Ca, ^^Zt, ^^^Sn and ^^^Sm) and where the diffuse surface plays an important role. In order 
to extend the analysis of ISGDR to these nuclei and also to exotic nuclei that are presently 
under intense investigations, one should adopt at first a more realistic assumption for the 
ground-state density distribution. In this case along with the density, other parameters of 
the nuclear Fermi liquid, e.g. the Lame coefficients acquire a radial dependence and the 
equations of motions must be solved numerically. 

In the present approach, the excitation of the ISGDR is not limited to the squeezing op- 
erator. It includes the entire momentum content in the operator ji{kr)Yi{f) and takes into 
account also its cm. correction via constraints on the density and displacement field fiuc- 
tuations. Moreover, since there is no mathematical or physical exception it considers along 
with the longitudinal solution also the transverse solution of the vector Helmholz equation. 
Conseqently the overtones of the ISGDR are mixtures of compressional and vortical velocity 
fields. For the second overtone, previously advocated to be of compressional nature, the 
nuclear Fermi-Fluid approach confirms very 

are pointing toward a cohabitation of compressional and vorticity vibrations in the ISGDR 
states up to 30 MeV. 

Naturally, a question arise : Has been gained so far any indication in experiment on the 
overtones with n > 3, located above the "high-lying" ISGDR state? For the time being this 
uestion cannot be answered because the latest data reported by the College Station group 
are providing very large uncertainties in the strength distribution for ^^'^Pb beyond the 
second peak, i.e. at excitation energy > Shw. 

A study of the electromagnetic multipole transitions was undertaken and concluded that 
the leading term in the Coulomb form-factor is singled-out only by the density fiuctuation 
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and can be related to the rms-charge radii. In turn the leading mutipole in the transverse 
electric form factor, the toroidal dipole moment, results soleley from the transverse part 
of the velocity field. In this respect the r.m.e. of the toroidal dipole transition provide a 
signature of vorticity through electromagnetic probes. 

The (e, e') are promising candidates for the exploration of the role of longitudinal and 
vortical currents of ISGDR since the separation of longitudinal (Coulomb) and electric trans- 
verse form-factor is feasible. However in such reactions it is difficult to avoid the excitation 
of the dominant electric dipole response, the IVGDR. It would be in this respect interesting 
to search for a macroscopical description that deals in a unified manner with the isoscalar 
and isovector electric dipole responses. 
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APPENDIX A: SCALAR AND VECTOR HELMHOLZ EQUATIONS 

The solution of the scalar HE (fTUI) reads 

V = J2(^xMkLr)Yxf.iO,4>) (Al) 

X/j, 

The solution of the vector HE splits into an electric(poloidal) and a magnetic (torsional) 
solution 

= E^aW^(^^^)^aa(^,0) (A2) 

Xfi 

^'"^ = E y=y {Syx-iVX+1 - 5a'a+iV^) jx'{kTr)Y^,A(^, 0) (A3) 

In the present study on electric resonances we are interested only in the poloidal solution. 
Using the properties of the spherical harmonic vectors we can derive the expression of the 
longitudinal and transverse displacements fields. Since 6sl results from the definition of the 
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scalar function 

V = V -Ss^V -Ssl (A4) 

we have that 



(A5) 



Similarly, from the definition of the vorticity 

1 



^poi = V X 5st (A6) 



we get 
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FIG. 5: Zeroth-order energy strengths distributions for p and J±. The energy cut is (£'*)cut 
MeV. 
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FIG. 6: First-order energy strengths distributions for p and J ±. The energy cut is (£'*)cut 
MeV. 
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